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Chapter  I 

Evidence  of  Correlation  in  Coupled  Systems 
l.A.  Introduction 

In  this  thesis,  two  separate  topics  in  queueing  theory  are 
discussed.  The  first  topic  concerns  the  quantitative  analysis 
necessary  to  support  intuitive  claims  that  correlation  exists  in  a 
system.  The  second  topic  investigates  the  simulation  of  certain 
dynamic  strategies  used  to  control  the  transmissions  in  a  time-slotted 
packet  radio  multiaccess  broadcast  channel  under  the  condition  that 
the  number  of  users,  N,  is  allowed  to  vary  slowly. 

The  dependence  due  to  access  protocols  intuitively  introduces 
correlation  into  a  system  of  nodes  in  a  packet  radio  network. 

However,  the  analysis  necessary  to  support  this  claim  has  been  found 
to  be  extremely  complex  and  tedious.  One  of  the  largest  problems  is 
simply  trying  to  produce  a  closed  form  expression  for  the  generating 
function  of  a  coupled  network.  This  thesis  shows  that  there  are 
methods  which  produce  results  that  support  the  claim  that  correlation 
does  indeed  exist  in  "intuitively"  correlated  systems. 

In  [6]  and  later  in  [7],  work  is  done  to  determine  average 
waiting  times  and  average  queue  lengths  for  two  node  systems.  A 
method  discussed  there  is  used  to  determine  conditional  and 
unconditional  expected  queue  lengths  for  symmetric  systems.  In  [1],  a 
method  is  introduced  to  determine  a  closed  form  expression  of  a 
generating  function  for  a  system  of  coupled  processors  by  reducing  the 
problem  to  a  Riemann-Hilbert  problem.  This  paper  gives  the  only 
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method  known  to  solve  for  the  generating  function  of  a  coupled  system. 
In  [8],  the  distribution  of  the  waiting  time  of  a  customer  is 
determined  for  a  shared  server  system.  A  general  background  for 
problems  exhibited  by  coupled  processors  is  found  in  [6],  and  a 
general  background  for  packet  radio  communications  where  much  of  this 
work  is  applicable  is  found  in  [5]. 

In  this  chapter,  the  discussion  of  correlation  begins  with 
describing  the  goals  of  this  work.  Models  are  presented  and  an 
efficient  method  is  used  to  determine  the  generating  functions  of  the 
models.  The  generating  functions  are  then  analyzed  to  produce  results 
consistent  with  the  goals  already  presented.  The  results  found 
concern  the  ordering  of  the  conditional  and  the  unconditional  moments 
of  the  systems,  the  ordering  of  the  poles  of  the  generating  functions 
given  for  some  systems,  and  the  phenomenon  that  hinders  the  search  for 
moments  of  higher  order.  An  ALOHA  model  is  then  presented  along  with 
a  set  of  goals.  In  this  model,  a  conjecture  is  supported  through 
simulation;  however,  a  valid  proof  of  this  conjecture  is  an  open 
problem.  This  is  followed  by  the  second  chapter  of  this  thesis, 
simulation  of  dynamic  strategies  used  to  control  transmission 
probabilities. 

l.A.l.  Goals  and  Models 

The  goal  of  this  chapter  simply  stated  is  to  produce 
quantitative  evidence  as  to  the  existence  of  correlation  in  a  system 
that  exhibits  correlation  intuitively.  At  the  basis  of  this  work  is 
the  idea  of  intuitively  correlated  versus  quantitatively  correlated. 
Although  one  system  may  seem  to  be  more  intuitively  correlated  than 


another,  a  measure  is  needed  to  make  such  a  claim  precise.  With 
quantitative  evidence,  such  claims  are  given  concrete  meaning. 

Two  such  quantitative  results  concern  the  generating  function 
determined  by  the  system.  The  first  is  a  relationship  between  the 
conditional  and  unconditional  moments  of  the  system.  If  the 
unconditional  expected  number  of  packets  in  a  station,  the  queue- 
server  combination  defined  later  for  the  models  in  this  paper,  is 
greater  or  less  than  the  conditional  expected  number  of  packets  in  a 
station,  given  that  the  other  station  is  empty,  correlation  exists 
between  the  two  stations. 

The  second  evidence  used  to  imply  correlation  concerns  the  poles 
of  the  generating  function  of  the  random  variable  determined  by  the 
total  number  of  packets  in  the  two  stations.  The  basis  for  the  argu¬ 
ment  that  poles  of  the  generating  function  are  a  possible  indicator  of 
correlation  lies  in  the  relationship  that  exists  between  the  pole 
placement  and  how  quickly  the  tail  of  the  distribution  function  goes 
to  zero.  It  is  seen  that  in  an  uncorrelated  system,  the  tail  goes  to 
zero  faster  than  in  a  correlated  system.  Intuitive  evidence  that  this 
is  true  is  found  in  the  argument  that  in  a  correlated  system  it  is 
expected  that  if  one  station  has  a  large  backlog,  then  the  other 
station  also  has  a  large  backlog.  This  results  in  the  random  variable 
determined  by  the  sum  of  the  packets  in  each  station  having  a  distri¬ 
bution  function  with  a  heavier  tail  than  the  tail  of  the  distribution 
function  of  an  uncorrelated  system.  A  heavier  tail  simply  means  that 
the  tail  goes  to  zero  more  slowly  than  a  comparable  tail.  This 


implies  that  an  ordering  of  the  poles  of  a  system's  generating  functions 
produces  evidence  of  correlation. 


Further  evidence  is  sought  involving  higher  order  moments  of  the 
generating  functions.  For  example,  if  the  first  joint  moment  of  the 
system  were  greater  than  the  product  of  the  first  moments  of  each 
station,  a  correlation  coefficient  could  be  determined  and  used  as 
evidence  of  correlation.  This  evidence,  however,  is  unattainable 
using  the  methods  presented.  The  problem  lies  in  the  fact  that  there 
are  more  unknowns  than  equations  available  to  be  solved  for.  In  fact, 
for  one  method  presented,  there  are  exactly  twice  as  many  unknowns  as 
there  are  equations  to  solve  them. 

An  emphasis  is  put  on  evidence  of  correlation  in  this  chapter; 
however,  the  discussion  of  higher  moments  is  felt  to  be  necessary 
since  it  is  directly  linked  to  the  methods  used  to  gain  results  for 
correlation. 

I.A.2.  Models  used  in  the  Analysis 

Three  different  systems  are  discussed  and  defined  here.  The 
first  is  a  system  which  will  be  called  the  splitting  system.  It  has 
customers  arriving  which  obey  a  Poisson  arrival  process.  Upon 
arriving,  the  customers  are  split  into  two  equal  packets,  and  each 
packet  is  then  queued  and  served  separately  by  two  stations.  The 
arrival  rate  is  A  and  the  service  times  are  exponentially  distributed 
with  rate  y.  If  each  station,  the  queue-server  combination,  is  looked 
at  separately,  the  stations  look  like  simple  M/M/1  queues  with  arrival 
rate  A  and  departure  rate  y.  This  system  exhibits  intuitive 
correlation  since  when  a  customer  arrives,  both  stations  receive  a 
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packet.  Therefore,  if  a  station  is  empty,  it  is  expected  that  the 
other  station  should  have  a  smaller  backlog  in  its  queue  than  if  the 
unconditional  backlog  is  looked  at.  Hence,  this  is  an  example  of  a 
correlated  system. 

The  next  system  discussed  is  similar  to  the  splitting  system 
in  that  the  stations  again  act  independently  upon  receiving  a  packet. 
The  difference  between  the  splitting  system  and  a  system  which  will  be 
called  the  independent  system  is  in  the  arrival  process.  In  the 
independent  system,  each  station  has  a  separate  independent  Poisson 
arrival  process  with  rate  X.  This  system  is  equivalent  to  two 
independent  M/M/1  queues.  The  results  from  this  uncorrelated  network 
are  compared  to  the  same  results  for  a  correlated  network. 

The  final  system  has  a  separate  independent  Poisson  arrival 
process  for  each  station  as  did  the  independent  system.  This  system, 
which  shall  be  called  the  shared-server  system,  uses  two  queues  and  a 
single  server.  The  server  has  exponentially  distributed  service 
times;  however,  the  rate  of  departures  is  2u.  The  server  alternates 
between  serving  the  queues  as  long  as  one  of  the  queues  is  not  empty. 
This  is  analogous  to  two  separate  servers  each  with  rate  u  for  each  of 
the  customer  types,  where  the  type  of  customer  is  determined  by  the 
queue  at  which  it  arrived,  if  both  customer  types  are  present  to  be 
be  served.  However,  if  one  of  the  queues  becomes  empty  and  hence 
one  of  the  customer  types  is  no  longer  present  in  the  system,  the 
single  server  will  put  its  entire  effort  to  serving  only  the  type 
of  customer  still  present  in  the  system  with  rate  2u.  Correlation 
exists  intuitively  in  this  system  since  if  only  one  customer  type  is 
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present  in  the  system,  the  server  will  work  with  rate  2u  rather  than 
with  rate  y  to  lower  the  backlog  in  the  nonempty  queue. 

The  shared-server  system,  although  it  is  an  intuitively 
correlated  system,  is  used  only  as  an  example  in  the  subsection  on 
determining  the  generating  function.  It  may  also  be  analyzed  using 
methods  similar  to  those  presented  in  this  paper  to  determine  evidence 
of  correlation.  However,  the  first  two  systems  presented,  the 
splitting  system  and  the  independent  system,  supply  ample  proof  that 
the  methods  used  to  provide  evidence  of  correlation  are  sufficient  in 
the  sense  that  the  goals  set  forth  are  satisfied. 

I.A.3.  Deriving  the  Generating  Functions 

This  section  uses  an  efficient  method  to  derive  the 
generating  functions  of  our  systems.  The  state  of  the  systems  is 
defined  to  be  the  number  of  customers,  or  packets,  in  station  i,  where 
i  =1,2  at  time,  t,  in  steady  state.  Consider  the  steady-state  joint 
generating  function  of  the  number  of  customers  in  each  of  the 
stations : 

n  (t)  n  (t) 

G(x,y)  =  lim  £[x  y  ] 

4 — K» 

where  n^(t)  and  n^t)  are  the  number  of  customers  in  stations  1  and  2, 
respectively,  and  the  Markov  chain  [n^(t),  n9(t)]  is  ergodic.  The 
determination  of  this  generating  function  is  the  goal  of  tm's 
subsection . 
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Letting 

|  =  Prob  (nx(t)  =  i,  n2(t)  =  j)  , 


D 


where  ir(i,j)  is  the  joint  steady-state  density  function  and  Q  is  the 
transition  rate  matrix  that  corresponds  to  the  system  under  study. 
These  conditions  are  necessary  for  an  ergodic  generating  function  to 
exist. 

The  transition  rate  matrix  may  in  some  cases  be  described  more 
simply  by  a  transition  rate  diagram  (see  Fig.  1).  This  is  true  for 
the  cases  given  by  the  systems  discussed  in  the  preceding  section. 

For  example,  the  transition  rate  diagram  given  in  Fig.  1.1  is  for  the 
splitting  system.  This  diagram  represents  customers  arriving  with 
rate  \  by  the  diagonal  transitions  and  customers  departing  with  rate  u 
by  the  horizonal  and  vertical  transitions.  Each  point  on  this  diagram 
represents  a  possible  state  in  which  the  process  may  be  found. 
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Fig.  1  Transition  Diagram  for  the  Splitting  System 

These  points  on  the  diagram  fall  into  four  categories,  points 
where  i>l,j>l,  points  where  i>l,j=0,  points  where  j>l,i=0,  and  the 
point  i=0,j=0.  These  categories  also  represent  four  different  parts 
of  the  generating  function  given  by  equation  (1.1).  These  parts  are 

[G(x ,y)  -  G(x,0)  -  G(0,y)  +  G(0,0)]  for  i>l,j>l, 

[G ( x , 0 )  -  G ( 0 , 0 ) ]  for  i>l,j=0, 

[2(0, y)  -  G(0,0)]  for  i=0,j>l,  and 

G(0,0)  for  i =0,j=0. 

Each  of  these  parts,  for  the  splitting  system,  has  a  different  arrival 
and  departure  process  assigned  to  it.  For  example,  the  point  i=0,j=0 
may  have  only  arrivals  whereas  the  point  i~l,jM  may  have  both 
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Under  the  condition-;  *yj.  '  ■  — ‘  t.  v>-  rr^nsi  tinn  diagram,  the 
allowable  transitions  for  the  fOur  categories  of  points  are  now 
incorporated  to  determine  the  generating  function  form  of  the 
equilibrium  relationship.  For  this  case  the  equilibrium  relationship 
is  0=  mQ.  For  the  splitting  system,  the  allowable  transition 
determined  by  an  arrival  introduces  the  term  '(xy-1)  since  an  arrival 
is  seen  as  an  arrival  to  both  stations.  A  departure  from  station  1  is 
represented  by  u(x~  -1),  and  a  departure  from  station  2  is  represented 
by  the  similar  result  u(y_1-l)-  Hence,  the  generating  function  form 
of  the  equilibrium  relationship  for  the  points  of  the  type  i>l,j>l  is 
simply 

[G(x,y)  -  G(x ,0)  -  G(0,y)  +  G(0,0)][\(xy-1)  +  u(x_1-l)  +  yfy"1-!)]. 

Using  this  technique,  the  entire  generating  function  form  of  the 
equilibrium  relationship  for  the  splitting  system  is  given  by 

0  =  [G(x ,y)  -  G(x ,0)  -  G(0,y)  +  G(0,0)][.\(xy-1)  +  u(x_1-l)  +  p(y_1-l)] 
+  [G(x,0)  -  G(0,0)][\(xy-1)  +  u(x_1-l)] 

+  CG(0»y)  -  G(0,0)][X(xy-l)  +  u(y_1-l)] 

+  G(0,0) [X(x.y-1 )]  . 

This  method  is  easily  applied  to  the  independent  system  where  the 
transitions  allowed  are  of  a  different  form  than  those  of  the 
splitting  system.  A  customer  arriving  at  this  system  introduces  terms 
of  the  form  \(x-l)  and  \(y-l)  depending  on  whether  the  customer 
arrives  at  station  1  or  station  2,  respectively.  The  generating 
function  form  of  the  equilibrium  relationship  then  takes  on  the  result 
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0  =  [G ( x ,y )  -  G( x ,0 )  -  G( 0 ,y )  +  G(0,0)][\(x-1 )  +  \(y-l)  +  u(x_1-l )  +  ^(y_1-l )] 
+  [G(x,0)  -  G(0,0)][\(x-1 )  +  X(y-l)  +  u(x  ^-1)] 

+  [G(0,y)  -  G(0,0)][\(x-1)  +  X(y-l)  +  y(y  *-l)] 

+  G(0,0 )[x(x-l )  +  X(y-l)]  . 

The  shared-server  system  also  has  different  transitions  allowed 
than  the  forementioned  systems.  For  the  points  ,j>l  and  the  point 
i=0,j=0,  the  transitions  allowed  are  the  same  as  in  the  independent 
system,  but  for  points  on  the  i=Q,j=0  axes  the  departure  rate  is 
increased  to  2y.  This  introduces  a  term  of  the  form  2u(x  ^-1)  or 
2y(y  *-l)  into  the  generating  function  form  of  the  equilibrium 
relationship.  This  result, 

0  =  [G(x,y )  -  G( x ,0 )  -  G(0,y)  +  G  (0,0)][X(x-l )  +  +  u(x_1-l  +  u(y_1-l)] 

+  [G(x,0)  -  G(0,0)][2u(x_1-1)  +  X ( x- 1 )  r  X(y-l)] 

+  [G(0,y)  -  G(0,0)][2y(y  *-l)  +  X(y-l)  +  X(x-l)] 

+  [G( 0,0 ) [\( x-1 )  +  X(y-1  )] 

is  consistent  with  the  previous  results. 

This  method  can  also  be  applied  to  discrete  time  processes  which 
have  an  equilibrium  form  of  m  =  op,  where  ~  is  previously  defined  and 
P  is  the  transition  probability  matrix.  However,  the  only  example  of 
this  type  of  system  is  given  in  section  l.C.  The  results  for  the 
so  lifting  system  and  the  independent  system  may  be  rewritten  as 


and 


G'U.y)  = 


(1.3) 


GI(x>y)  =  GI(x,0)(y~1-l)  +  GI(0,y)(x~1-l) 

(x-1- 1 )  +  (y_1-l)  +  p(x-1)  +  p(y-l) 

C  T 

where  o  =  X/u  and  G  (x,y)  and  G  (x,y)  are  the  notations  used  to  denote 
the  generating  function  of  the  splitting  system  and  the  independent 
system,  respectively.  These  two  results  are  now  used  to  show  that  the 
goals  set  at  the  beginning  of  this  chapter  may  be  met. 


l.B.  Indicators  of  Correlation 


In  this  section,  the  ordering  of  the  conditional  and  of  the 
unconditional  moments  of  first  order  of  the  systems  is  investigated. 
Following  this  discussion,  the  other  evidence  of  correlation  is  then 
presented  using  an  argument  which  depends  on  the  ordering  of  the  poles 
of  a  generating  function.  Higher  moments  of  the  systems  are 
investigated  as  the  final  topic  of  this  section. 


I.B.l.  Ordering  of  the  First  Order  Conditional  and 
Unconditional  Moments 

Both  the  independent  and  splitting  systems  are  analyzed  to 
determine  the  ordering  of  their  higher  order  moments.  Equations  (1.2) 
and  (1.3)  are  multiplied  by  xy  in  the  numerator  and  in  the  denominator 
to  simplify  the  analysis.  Some  simple  properties  of  the  generating 
functions  are  new  stated  which  will  be  used  later  in  the  analysis. 
Since  the  systems  stuaied  are  symmetric, 


P'-ccerty  1 
Prcp0^*tv  2 


G(x,y)  =  G(y ,x )  ,  and 


a  ,  ( x ,  v  ;  =  G  , 

k.on  x  '  '  ,iTuk 

L  w  1  £ 


1 


(y,x)  ,  where 


12 


,k 


jn 


V,- (x>y)  riG<x'y) 

12  ax  ay 


are  both  properties  of  the  generating  function.  Another  important 
property  i s : 


Property  3  Tim  G  ,  (x,l)  =  lim  G  ,  (x,x) 

x-1  lk  x-1  1K 


This  is  a  consequence  of  the  analyticity  of  the  generating  function. 
The  final  property  presented  follows  directly  from  calculus: 

Property  4  [G(x,x)]  =  G1  (x,x)  +  G2  (x,x) 

Using  these  properties,  the  conditional  and  unconditional  first 
order  moments  are  now  determined  for  the  systems  under  investigation. 
Substituting  in  y=l  for  the  altered  (1.2)  as  described  above,  and 
using  property  1,  the  result  that 

gS(x,D  =  "T-XO^  (1-4) 

is  attained.  Similarly,  for  the  independent  system,  the  result  that 

g'ix.!)-  (!-5) 

is  found.  Using  the  common  normalizing  condition,  G(i,l)  =  1,  the 
result  that 

GS( 3 , 1 )  =  G1 (0,1 )  =  (1-c) 


1»  i  . 


i 


provides  the  probability  of  one  of  the  stations  being  empty. 
Expanding  the  results  from  (1.4)  and  (1.5),  the  equation  that 
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I 

I 


G^x.l)  =  GS( x ,1 )  = 
dxK  dxk 


(l-px)k+1 


(1.6) 


is  produced.  When  k  is  allowed  to  equal  1  and  x=l ,  the  unconditional 
expected  number  in  a  station  is  given  by  (1.6). 

To  derive  the  conditional  expected  number  in  a  station  given  that 
the  other  station  is  empty,  the  derivative  of  G(x,y)  given  by  (1.2) 
and  (1.3)  for  each  system  is  taken  with  respect  to  x  (or  y).  This 
result  is  followed  by  setting  y=x  which  produces  the  following  for  the 
splitting  system: 

C  G^(x,0)  GS(x,0)p( 2x+l) 

G i  (x,x)  =  - - -  +  - 5-  (1.7) 

(2-xp(l+x))  (2-xo(l+x); 


and  for  the  independent  system 

r  G, !(x,0)  2eGI(x,0) 

G, 1  (x,x  )  *  — -  +  - 5-  '  (1.3) 

(2-xo)  (2-xoT 

Allowing  x-1 ,  G,(1,0),  which  is  equal  to  G?(0,1)  by  property  2,  is 
attained  for  both  systems: 

G^d.O)  =  c/2  ,  and 
G,  1  ( 1 ,0)  =  p  . 


E 


These  results  are  used  to  gain  tne  expected  number  in  one  station 
given  that  the  other  is  empty,  since  this  value  is  equal  to  6^1 ,0)/G(l  ,0) 
Hence,  the  conditional  first  order  moment  for  one  station,  given  that  the 
other  station  is  empty,  is  p/2(i-p)  for  the  splitting  system  and  o/(l-o). 
for  the  independent  system.  These  results  provide  the  first 
quantitative  evidence  that  the  splitting  system  is  correlated  and  that 
the  independent  system  is  uncorrelated .  Hence,  the  conditional  moment 
is  less  than  the  unconditional  moment  in  the  splitting  system  which 
corresponds  to  the  intuitive  claim  that  if  one  station  is  empty,  the 
other  station  tends  to  have  a  smaller  backlog  than  when  the  oacklog  is 
not  conditioned  upon  having  one  station  empty.  In  the  independent 
system,  both  moments  are  equal,  implying  that  conditioning  upon  an 
empty  station  does  not  affect  the  backlog  in  the  other  station,  and 
that  the  system  is  uncorrelated. 

This  technique  provides  only  one  piece  of  evidence  as  to  the 
existence  of  correlation.  Further  insight  into  the  techniques 
available  to  produce  evidence  of  correlation  is  necessary  and  given  in 
the  next  subsection. 

I.B.2.  Ordering  of  Poles 

This  subsection  discusses  the  ordering  of  the  poles  of  the 
generating  function  for  the  random  variable  defined  as  the  sum  of  the 
number  of  customers  in  stations  1  and  2  for  the  splitting  system  and 
rcr  the  independent  system.  A  third  system  is  also  defined  to  provide 
another  pole  to  compare  to  the  poles  already  found. 

This  system  is  an  intuitively  totally  correlated  system  defined 
by  two  stations  using  an  arrival  process  similar  to  that  of  the 
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splitting  system.  However,  the  customers  must  also  depart  at  the  same 
time  unlike  either  the  splitting  or  the  independent  systems.  This 
simply  says  that  if  a  packet  arrives  at  the  network,  it  appears  as  a 
customer  arriving  at  each  station.  If  a  customer  departs  from  one 
station,  there  is  always  a  departure  from  the  other  station.  The 
arrival  process  is  Poisson  with  rate  X, and  the  departure  process  has 
an  exponentially  distributed  service  time  with  rate  y.  Using  the 
technique  described  in  subsection  I.A.3.,  the  joint  steady-state 
generating  function  is  determined  to  satisfy 

G^(x,x)  =  where  p  =  X/y  (1-9) 

1-px^ 

for  y=x.  This  is  sufficient  since  the  generating  function  of  the 
random  variable  defined  by  $e  sum  of  the  number  of  customers  in  each 
station  is  equivalent  to  setting  y=x  in  the  joint  steady-state 
generating  function.  Another  way  of  looking  at  this  system  is  to  view 
it  as  just  a  M/M/1  queue  which  allows  only  two  customers  to  arrive  or 
to  depart  instead  of  just  one  customer  arriving  or  departing.  This 
system  can  be  interpreted  s  totally  correlated  since  each  station 
contains  the  same  number  of  packets. 

The  tie  between  the  poles  of  the  generating  function  of  the 
random  variable. which  is  defined  by  the  sum  of  the  number  of  customers 
in  each  station  denoted  by  Z  and  correlation, lies  in  the  relationship 
between  the  poles  and  the  tail  of  this  random  variable.  It  is  known 
that  for  station  i,  i  equalling  1  or  2,  each  of  the  three  systems 
discussed  has  the  same  marginal  distribution.  Therefore,  a  heavier 


tail  for  the  random  variable  given  oy  the  sum  of  the  number  of 
of  packets  at  each  station  indicates  that  more  correlation  exists 
in  that  system.  So,  by  determining  if  the  tail  for  one  system  is 

heavier  than  that  of  another,  evidence  is  obtained  supporting  the 
claim  that  more  correlation  exists  in  one  system  than  in  another. 

One  way  to  measure  how  heavy  a  tail  is  is  to  look  at  the  least 
magnitude  pole  of  the  generating  function  of  Z.  For  example,  for  the 
geometric  distribution  with  parameter  3,  the  pole  of  the  generating 
function,  given  by  1/3,  is  inversely  proportional  to  the  magnitude  of 
3.  This  results  in  a  heavier  tail  for  a  larger  3,  or  equivalently  for 
a  smaller  least  magnitude  pole.  In  general,  roughly  speaking,  if  for 
one  system,  the  pole  of  least  magnitude  is  found  to  be  less  than  that 
for  another  system,  the  first  system  has  a  tail  that  goes  to  zero 
faster  than  the  second  system's  tail. 

Using  equations  (1.2)  and  (1.3),  where  y  is  set  equal  to  x,  and 
(1.9),  the  poles  of  the  generating  function  of  the  random  variable,  Z, 
are  found  to  be: 

x!  =  l/o 

x'^  =  1/v’p  ,  and 


where  x“,  j=I,  TC,  or  S,  corresponds  to  the  poles  of  the  independent, 
totally  correlated,  and  splitting  systems,  respectively.  The  poles  of 
totally  correlated  and  independent  systems  are  easily  determined. 
However,  the  poles  of  the  splitting  system  are  also  affected  by  the 

C 

location  of  any  possible  poles  of  the  unknown  G  ( x , 0 ) .  And,  for  x  real. 
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GS(0,x)  <  GS(l,x) 

s 

and  the  term  G  (l,x)  is  known  to  have  a  first  pole  at  1/p,  so,  since 

c  c  c 

G  (0,x)  can  not  go  to  infinity  prior  to  G  (l,x),  G  (0,x)<+®  for  0<x<l. 

C 

So, any  pole  of  G  (0,x)  does  not  interfere  with  the  placement  of  the 
first  major  pole  for  the  splitting  system,  and  the  ordering  of  the 
poles  follows  as 

xTC  <  xS<  x1  for  pe(0,l)  . 

So,  by  using  the  tail  argument,  this  ordering  implies  that  the 
splitting  system  is  more  correlated  than  the  independent  system,  and 
less  correlated  than  the  totally  correlated  system. 

This  technique  provides  the  second  piece  of  evidence  used  to 
show  the  existence  of  correlation  in  a  system.  Further  evidence  is 
sought  in  the  next  section,  but  no  definitive  results  are  produced. 

I.B.3.  Higher  Order  Moments 

In  this  subsection,  relationships  between  the  higher  moments 
are  given  and  a  general  method  to  determine  these  relationships  is 
also  discussed.  Although  the  ultimate  goal  is  to  solve  for  these 
moments  explicitly,  this  goal  is  unattainable  using  the  methods  given 
here  at  this  time.  However,  it  is  of  interest  that  using  a  method 
presented  here,  only  one  piece  of  information  is  gained  for  every  two 
iterations  of  the  method.  Two  methods  are  discussed  in  this 
subsection.  The  first  is  an  extension  of  the  work  done  in  subsection 
l.B.l.  while  the  second  uses  the  analyticity  of  the  denominator  of  the 


generating  function  around  the  point  x=l,y=l  on  the  curve  defined  by 
the  zeroes  of  the  denominator.  When  the  sets  of  equations  found  using 
these  two  methods  are  combined,  there  are  still  not  enough  equations 
to  solve  for  any  of  the  higher  order  unknowns.  Both  methods  are  applied 
only  to  the  splitting  system  here.  However,  other  systems  may  also 
be  analyzed  using  these  methods. 

The  first  method  discussed  is  a  direct  extension  of  the  analysis 

done  in  subsection  l.B.l.  The  same  properties  hold,  and  the  results 

found  there  are  necessary  here  to  continue  this  analysis.  Using  the 

5 

result  already  found  for  (x,y),  the  derivative  with  respect  to  y  is 

c 

taken.  Again,set  y=x  and  allow  x+1.  The  result  for  G^  (1,1)  is  the 
following : 

GS2  (1,0)  (p2+p) 

G  S(i,i)=-J - + - j  .  (1.10) 

2 ( 1- p )  2(l-p)2 

This  equation  has  two  unknowns,  so  it  is  not  solvable  for  either 
unknown.  Further  derivates  only  increase  the  number  of  unknowns  by 
two  for  each  new  equation,  namely,  two  new  joint  moments  of  the  form 

I1-1)  a"d  GjkL  f1-0'  • 

This  does  produce  a  set  of  (k+m-1)  equations  in  2(k+m-l)  unknowns  which, 
upon  the  introduction  of  (k+m-1),  more  independent  equations  in  the 
same  unknowns  are  solvable. 
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The  second  method  to  determine  higher  order  moment  equations 
revolves  upon  analyzing  the  zeroes  of  the  numerator  of  the  generating 
function  in  a  neighborhood  around  x=l,y=l  on  the  curve  determined  by 
the  zeroes  of  the  denominator  of  the  generating  function.  The 
splitting  system  is  used  as  an  example  for  this  method,  and  equation 
(1.2)  is  altered  by  multiplying  the  numerator  and  the  denominator  by 
xy  to  produce 


GS(x,y) 


H(x)x(l-,y)  +  H(y)y(l-x) 
x(l-y)  +  y(l-x)  +  pxy (xy-1 ) 


(1.11) 


s  s 

where  H(x)  =  G  (x,0)  and  H(y)  =  G  (G,y).  This  is  a  much  simpler  form 
of  the  generating  function,  and  it  is  much  easier  to  analyze  in  this 
form.  Now  the  curve  defined  by  setting  the  denominator  to  zero  is 

y(l-x)  +  x(l-y)  +  pxy (xy-1 )  =  0  .  (1.12) 


The  point  x=l,y=l  is  a  solution  to  this  equation,  hence  it  lies  on  the 
curve  determined  by  equation  (1.12).  Substituting  x=l,y=l  into  the 
numerator  of  equation  (1.11),  it  is  found  that  this  point  is  also  a 
solution  to  the  equation, 

xH(x ) ( 1-y )  +  yH(y)(l-x)  =  0  .  (1.13) 

Therefore,  the  point  x=l  ,y=l  results  in  an  indeterminate  form  for 

s  s 

G  (1,1).  It  is  known  that  G  (x,y)  is  a  generating  function,  so  it  is 
analytic  in  a  ball  about  the  origin  which  contains  the  point  x=l,y=i. 
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Using  this  knowledge,  a  solution  is  found  for  equation  (1.13)  for  a 
function  y  =  i(x).  This  result, 

H(*(x))  =  (llj/xT  H(x)  U-14) 

and  the  result  that  is  found  when  equation  (1.12)  is  solved  for  y=i(x), 
4>(X)  =  -  - JH4-L  (1.15) 


( l-2x-ox )  -  ((l-2x-ox)2  -  4cx3) 
2ox2 


produces  a  foundation  for  this  method.  By  solving  for  f(x)  in  equation 
(1.15)  and  then  substituting  this  result  into  equation  (1.14),  it  appears 
it  is  possible  to  solve  for  H(l)  and  for  all  the  derivatives  of  H(x) 
with  respect  to  x  evaluated  at  H(l).  Before  proceeding  with  this  method, 
a  simpler  format  is  presented  to  simplify  the  necessary  calculations. 
Using  the  generating  function  defined  as, 

where  N(x,y)  is  equivalent  to  equation  (1.13)  and  D(x,y)  is  equivalent 
to  equation  (1.12)  and  using  the  fact  that  G  (x,y)  is  analytic,  so  when 
D(x,y)=0,  the  equation,  N(x,y)=0,  is  also  true,  the  method  produces 
results  for  N(x,y).  This  is  possible  since  using  the  ;(x)  defined  in 
equation  (1.15)  and  using  the  result  that, 

N ( x ,  :•  ( x ) )  =  0  and  D(x,  ;(:<))=  0 

as  well  as  the  result  that  all  the  derivatives  of  N(  x,  :(;*))  and  ^ 

D(x,  :(:<))  are  equal  to  zero.  N'(x,  :(;<))  and  all  the  derivatives  of 
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m 

fk  tV  N(x’y)  =  N  k  m  • 

3x  3y  x  y 


„k 

7T  zs  0(x’y)  =  V,.  ■ and 

ox  o y  x  y 


„k  k  times 

77  ?(*)  =  ■*(*) . 


=  ?’  k(x) 


it  is  possible  to  produce  tables  of  solutions  for  these  values.  For 

example,  Table  (1.1)  is  a  list  of  all  the  possible  D  .  m  values 

k  x  y 

for  tr.  splitting  system  with  x=l,y=l.  The  *>'  (1)  values  are 

k 

determined  by  taking  derivatives  of  D(x,p(x))  and  by  solving  for  (x) 
with  x=l.  For  example,  taking  the  first  derivative  of  D(x,s(x)) 
results  in 


D  +  (x)D  =  0  . 

x  y 


Table  (1.1)  Table  of  D  k  m  values  for  the  point  x,y=l 

x  y 

of  the  splitting  system 


0(1,1)  =  1 

0x(i,l)  =  Dy  (1,1)  =o-l 
D  ,(1,1)  =  D  ,(1,1)  =  2o 


VM)  =3,-2 


n  ( 1  M  =  n  m  i  \  =  1 


J  2 v  *  ’  '  ?  1 

y:<"  yx" 


0  9  -,(1,1)  =  4  o 


0  k  _  =  3  for  (r?=0,  k>2),  (.<=0,  r;>2\ 
■<' /’  (~.=i  ,  k>2  ; ,  (k  =  l,  •:>:), 


and  (m>2,  k>2; 


Solving  for  *'(1),  the  result  is  j * ( 1 ) = - 1 .  Continuing  this  procedure 
a  table  of  5'  (1)  is  computed  and  compiled  in  Table  (1.2)  for 
ke{0, 1,2,3}. 

Table  (1.2)  Table  of  <fi,k(l)  values  for  ke{0,l,2,3} 


b(l)  =  1 
<D'(1)  =  1 

r2d)  - 


4  -  2: 


b,J(D  « 


(l-o)‘ 


A  table  for  N  k  where  x=l,y=l  is  presented  in  Table  (1.3). 
x  y 

Table(l.3)  Table  of  N  ,  (1,1)  for  the  splitting  system 

x  y 


Nx(l,l)  -  Ny (1,1)  =  -H( 1 ) 

N  2(1 ,1)  =  N  2(1 ,1 )  =  0 
x  y 


N  (1,1)  =  -2 (H( 1 )  +  H'(l)) 

AJ 

N  2(1 ,1)  =  N  2(1,1)  =  -H ' ( 1 )  -  H'2(l) 


N  k  irt  ( 1  ’  1 )  =0 


for  k  ,j>2 


N  k  (1,1)  =  N  k  (1,1)  =  -kH 1 k" 1 ( 1 )  -  H ' k ( 1 ) 
x  y  y  x 


Lisina  a  similar  method  as  was  used  to  solve  for  :(:<)  and  for  the 


derivatives  of  :(;<)  by  the  solution  ta  the  equations  using  tne  D(x,y 
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terms,  it  is  possible  to  produce  an  equation  in  H(x)  and  in  the 
derivatives  of  H(x)  by  using  the  N(x,y)  terms  in  place  of  the  D(x,y) 
terms.  For  example,  the  first  derivative  of  N(x,y),  where  y=  $(x), 
produces , 


Nx  +  Ny  to’  (1)  =  0  . 


(1.16) 


Using  the  values  found  in  Tables  (1.2)  and  (1.3),  it  is  seen  that 

N  =-N  which  is  just  an  identity.  Consequently,  from  this  point  on 
x  y 

all  the  N(x,y)  terms  are  evaluated  at  x=l ,  y=l,  so  the  (1,1)  has  been 
dropped.  The  next  derivative  of  equation  (1.16)  produces 


N  +  2  N  <j> '  ( 1 )  +  N  (4> '  ( 1 ) ) 2  +  N  4> '  2  ( 1 )  =  0 
xx  xy  v  v  ‘  yyvv  v  "  y 


(1.17) 


Solving  this  for  H 1 ( 1 )  using  Tables  (1.2)  and  (1.3)  results  in 
H 1  ( 1 )  =  0/2  . 


This  is  the  same  result  found  in  subsection  l.B.l.  for  G^(1,0).  The 
next  derivative  of  equation  (1.17)  produces  another  identity  as  was 
the  case  for  equation  (1.16).  The  next  derivative  following  this  one 
results  in  the  equation 

N.+4N  .  V(l)  +  6  N  9  ?(i'(l))2  +  6  N  9?'2(1)+4N  ,(:'(!): 


+  12  N  ( ; '  (1 ) )  ( :/ 2(  1 ) )  +  3  N  2(:’2(1))2+  4  N  2  ( : '  ( 1 ) )  ( : '  3  ( 1 ) ) 
^  v  x1- 


+  Ar{:)f  +  5  n  ,(:'(i))2(: |2(D)  +  4  :j 

y  yJ  y 


•  '-('ll 
\  1  / 


■Ny  :  ‘ 4  ( : 


(1.13) 


wnich  simplifies  In 


+  12  N  2(«y(l))U,2(l))  +  4  N  ^3(1)  +  N  i'4(l)  =  o  (1.19 
xy  y  y 


upon  the  substitution  of  the  terns  known  to  be  zero  from  Table  (1.3). 
Now,  using 


y4(1)  =  24(8-llo  +  6o2  -  p3) 

(1-p)  3 


and  Table  (1.3),  it  is  found  that  equation  (1.19)  reduces  to  the 
equation , 


AH  '  3 ( 1 )  +  BH,2(1)  +  C  =  0 


(1.20 


where  A  =  3, 


B  =  24  + 


and 


c  .  (  4(1)  ,24^,114^ 

(1--')^  (1-.) 


Upon  taking  further  derivatives,  each  even  numbered  derivative  seems 
to  produce  an  identity  and  each  odd  numbered  derivative  seems  to 
result  in  another  equation  in  two  more  unknowns.  This  trend  is 
expected  to  continue  for  all  derivatives,  but  it  has  been  veri- 
Tied  only  up  to  the  eighth  derivative. 

These  two  methods  produce  two  separate  sets  of  eouations  that 
independently  are  not  sufficient  to  solve  for  any  of  the  higher  order 
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moments.  Although  it  is  strai ghtforward  to  solve  for  the  first  order 
moment  due  to  the  symmetry  that  exists  in  these  systems,  the  necessary 
independent  equations  needed  to  solve  for  the  higher  order  moments 
are  lacking.  Together,  the  two  sets  of  equations  in  the  best  case 
produce  two  equations  in  three  unknowns.  In  the  general  case  when 
the  first  method  has  (k+m-1)  equations  in  2 ( k+m- 1 )  unknowns,  the 
greatest  number  of  independent  equations  the  second  method  supplies  is 
(k+m-1 )/2  which  is  far  short  of  the  necessary  (k+m-1)  equations  needed 
to  solve  for  any  of  the  unknowns.  For  this  reason,  the  solution  to 
higher  order  moments  is  still  an  open  problem. 


l.C.  Evidence  of  Correlation  in  a  Three-Node  Slotted  ALOHA  Network 

This  section  uses  simulation  to  provide  evidence  that  correlation 
exists  in  a  three-node  slotted  ALOHA  network.  As  a  direct  extension 
of  the  work  done  in  [9]  for  a  two-queue  network,  if  the  inequality 
shown  to  be  true  by  simulation  is  always  true,  the  thrcugnput  for  a 
three-queue  buffered  ALOHA  network  can  also  be  determined  [2]. 


l.C.l.  Model  and  Goals 

The  model  used  is  a  three-node  slotted  ALOHA  broadcast 
system.  Stations  one  and  two  have  transmission  prcbabi 1 i ties  f^  and  f„ 
when  they  nave  a  oacket  to  transmit,  respectively,  and  mean  arrival 
rates  a.  and  a ? ,  respectively .  The  third  station  has  an  infinite 
backlog  and  transmits  with  orobability  p  in  each  slot.  The  state  of 


ied  o / 


on  one  an( 


nur.cer  oi 


in  station  two  at  the  end  of  3 


the  methods  discussed  in  subsection  I.A.2.,  a  generating  function  form 
of  the  equilibrium  equation,  m=  -tP,  where  -r  is  the  joint  steady-state 
density  function  and  P  is  the  transition  probability  matrix,  is 
found  to  be 

G(x,y)  =  A(x)A(y)  jjG(x,y  )  -  G(x,0)  -  G(0,y)  +  G(0,0 )  ]B(  x  ,y  ) 

+  [G(x ,0 )  -  G ( 0 ,0 ) ]C( x )  +  [G(0,y )  -  G(0,0)]C(y) 

+  G  ( 0 , 0 )] 

where 

B(x,y)  =  (1  +  pfjf^x'1-!)  +  ‘pf2f1(y“1-l)) 

C(x)  =  (1  +  pf 1 (x-1-l ) )  , 

C(y )  =  (1  +  pf2(y“1-l))  ,  and 

p  =  (1-p),  ft  =  (1-fj),  and  f?(l-f2) 

The  A(x)  and  A(y)  terms  represent  the  generating  function  of  the 
arrival  process.  In  the  model  discussed  here,  they  are  represented  by 
a  Bernoulli  arrival  process  with  rates  a^  and  a^  so 

A(x )  =  1  +  a^ (x-1 ) ,  and 
A  (y )  =  1  +  a2(x-l) 

are  the  generating  functions  of  the  arrival  process. 

The  goal  of  this  work  is  to  show  by  simulation  that  the 
probability  of  being  in  the  state  where  stations  one  and  two  a>le  empt;. 
is  less  than  the  product  of  the  probability  that  station  one  is  emoty 
and  tne  orobaoility  that  station  two  is  empty.  Written  another  way , 
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G(0,0)  G(0,1)G(1 ,0)  (1.21) 

The  simulation  determines  estimates  of  6(0,0),  G(0,1),  and  G(1,0). 


I.C.2.  Simulation  and  Results 

The  simulation  uses  a  set  of  conditions  to  determine  the  a^, 
a2»  f2»  and  P  necessary  to  determine  the  unknowns  in  equation 

(1.21).  The  program  used  is  written  in  the  Fortran  language  and  is  in 
no  way  considered  to  be  efficient.  The  program  simulates  arrivals  and 
departures  in  an  ALOHA  system.  The  structure  is  as  follows.  It  uses 
a  main  program  to  generate  random  variables  and  to  call  the  first  sub¬ 
routine.  The  random  numbers  chosen  between  0  and  1  are  used  to  deter¬ 
mine  if  an  arrival  or  a  departure  took  place  later  in  the  program. 

The  first  subroutine  sets  the  parameters  mentioned  earlier  such  that 
they  satisfy  the  following  conditions. 


p  =  ( .  1 ) i  -  0.05 
f1  =  (.l)j  -  0.05 

f2  =  ( . 1 ) k  -  0.05 

a^  =  (. l)m  -  0.05 

a2  =  ( . 1 )n  +  0.05 

ai  <  f1(l-f2)(l-o) 

a 9  ^  *2'  ^~r  1 


for  i=l,10 
for  j=l  ,10 

for  k=l  ,10 

for  m=l  ,10 

for  n  =  l  ,10  ,  and 


The  last  two  conditions  limit  the  region  that  the  simulation  is 
computed  for  by  restricting  the  parameter  to  fall  within  a  strongly 
staple  region. 


The  second  subroutine  uses  the  passed  parameters  to  stimulate  an 
ALOHA  system.  A  fifty-thousand  slot  run  is  made  to  compute  the  — 

estimated  number  of  times  the  system  is  in  the  states  described 
earlier.  An  arrival  occurs  in  a  station  if  the  parameter  value  is 
larger  than  an  independent  random  number  between  0  and  1  produced  in 
the  main  program.  A  transmission  occurs  if  the  station  is  not  empty, 
no  other  station  is  transmi tti ng,  and  the  parameter  value  for  that 
station  determined  in  the  first  subroutine  is  larger  than  the 
independent  random  number  between  0  and  1  tabulated  in  the  main 
program.  The  states  of  both  stations  are  kept  track  of,  and  if  either 

one,  the  other,  or  both  stations,  are  empty,  a  flag  is  set  corres-  ~ 

ponding  to  the  condition  present  during  that  slot.  The  number  of 
times  a  certain  flag  is  set  is  kept  track  of  and  after  the  50,000 

slot  run,  the  number  of  times  a  specific  flag  is  set  is  divided  by  .j 

50,000.  These  quantities  are  estimates  of  G(0,0),  G(1,0),  and  G(0,1). 

A  printout  of  the  estimates  of  G(0,0)  and  the  product  G(0,1)G(1,0)  is 
used  to  verify  equation  (1.21)  for  the  parameter  values  used.  ^ 

Over  the  total  range  of  the  parameters  used,  equation  (1.21)  held 
for  the  estimates.  This  does  not  say  that  equation  (1.21)  is  true  for 
all  parameter  values,  but  it  does  produce  hope  that  a  proof  of  equation  *■» 

(1.21)  is  possible. 
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Chapter  II 


Verification  that  Predictions  are  Valid  for 
Control  Strategies  Using  Simulations 

2. A.  Introduction 

In  this  chapter,  a  time- slotted  ALOHA  broadcast  channel  is 
simulated  using  two  different  types  of  channel  feedback  to  determine 
transmission  probabilities.  The  simulations  done  are  an  extension  of 
the  wor\  done  in  [3]  where  strategies  based  on  stochastic  approxima¬ 
tions  are  presented  to  control  transmission  probabilities  of  stations 
in  a  collision  access  communication  channel.  The  main  change 
incorporated  into  the  simulation  is  that  the  number  of  users  is 
allowed  to  "slowly"  vary. 

In  this  section,  yoals  and  a  review  of  the  work  done  in  [3]  are 
presented.  We  describe  in  section  B  the  simulations  of  the  fixed  user 
model  with  a  nonzero  e,  a  parameter  that  in  [3]  is  allowed  to  tend 
toward  zero.  The  data  from  these  simulations  are  then  compared  to  the 
values  given  in  [3].  In  section  C,  the  variable  user  models  are 
presented.  Simulations  of  these  models  are  also  done  and  the  data 
from  these  simulations  are  reviewed. 

2.A.I.  Goals 

The  motivation  of  this  work  is  to  simulate  the  0-1-e 
feedback  and  the  acknowledgment  based  feedback  strategies  under  the 
conditions  that  the  parameter  e  is  strictly  greater  than  zero  and 
that  the  number  of  users  is  allowed  to  slowly  vary  in  order  to 
verify  that  the  theory  presented  in  [3]  is  capable  of  predicting 
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simulation  results.  We  want  to  see  that  under  conditions  where  the 
theory  cannot  be  applied  directly  to  produce  exact  analytic  results, 
as  in  both  of  the  forementioned  cases,  it  is  still  possible  to  predict 
the  behavior  of  the  control  strategies. 

2. A. 2.  Background  Information 

This  subsection  reviews  work  first  presented  in  [3].  Only 
work  relevant  to  the  simulations  done  in  this  paper  is  presented.  The 
models  used  throughout  this  paper  consider  N  stations  sharing  a  com¬ 
munication  channel,  where  N  is  not  assumed  to  be  fixed  unless  stated 
so.  The  channel  can  be  accessed  only  during  specific  time  slots  which 
result  in  a  station  j  of  the  N  possible  stations  transmitting  a  packet 
during  time- slot  k  with  probability  f.(k).  If  station  j  transmits  in 

J 

slot  k,  U. (k)  equals  1  and  U . ( k )  equals  zero  otherwise.  Station  j 

J  J  ' 

receives  channel  feedback  at  the  end  of  a  slot.  This  feedback  is 
represented  by  Y.(k).  It  is  assumed  that  given  £(k),  the  random 

J 

vector  given  by  (U(k),  _Y(k))  is  conditionally  independent  of  the  past 
and  that  the  conditional  distribution  of  (JJ(k),  _Y(k))  given  jf(k)  is 
not  dependent  on  k.  The  transmission  control  strategies  considered 
take  on  the  form 

fjU+l)  =  [f  j(k)  +  e  G j ( f j ( k ) ,  Yj (k) ,  UJ(k))]^J'  (2.1) 

where  A  =  (A^ ,  ...  ,  A^ )  is  a  vector  with  positive  coordinates,  [x]^ 
denotes  the  projection,  min  { b ,  max  { a , x } } ,  of  x  onto  the  interval 
[a,b],  and  the  G.'s  are  bounded  continuous  functions. 

By  the  given  assumptions,  there  is  an  N-vector  valued  function  G 


defined  by 


(2.2) 


GjCf)  =  E  CSj(fj(k),  Yj(k).  Uj(k))|f(k)  =  f]. 

By  allowing  a  time  transformation  where  ke  =  t  for  an  integer  k,  it  is 
possible  to  define  a  continuous- time  random  process  fe  by  setting 
fe(t)  =  f(k)  where  it  is  also  required  that  the  sample  paths  of  fe  be 
affine  over  intervals  of  the  form  [ke,  (k+l)e].  This  is  equivalent  to 
expanding  time  using  the  parameter  t  while  e  tends  to  zero.  Then 
using  the  theory  of  stochastic  approximation,  under  general 
conditions,  the  random  process  f£  converges  weakly  as  e  tends  to  zero 
(essentially  all  that  is  needed  is  that  the  conditional  distribution 
of  (]J(k),  _Y(k))  given  f(k)  should  depend  smoothly  on  f_,  see  [4, 
sections  5.3,  5.4,  and  4.4])  to  the  solution  of  the  following  ordinary 
differential  equation  (Q.D.E.)  with  boundary: 

f  =  6(f)  +  v.  +  u  •,  1(0)  given 

0  <  f(t)  <  A  ,  v(t)  <  0,  u(t)  >  0  (2.3) 

vi(t)(A.  -  f.(t))  =  0  ,  u.(t)  ^(t)  =  0 

where  _u  and  v^  are  auxiliary  functions  to  be  determined  together  with 
!•  This  convergence  result  is  akin  to  the  law  of  large  numbers.  A 
second  convergence  result,  akin  to  the  central  limit  theorem,  pertains 
to  the  limiting  distribution  as  e  tends  to  zero  of  the  random  process 
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where  f  is  a  locally  asymptotically  stable  equilibrium  point  of  the  ODE 

{ 2.3 ).  To  simplify  the  analysis,  it  is  assumed  that  0  <  f .  <  A.,  for 

J  J 

each  j.  Then  under  general  conditions  (see  [4,  sections  5.6  and 
4.7]),  if  f£(0)  converges  weakly  (as  e  tends  to  zero)  to  a  constant 

if 

x(0),  or  if  f  is  globally  asymptotically  stable  and  on  initial 
transient  period  is  ignored,  f£  converges  weakly  to  a  Gaussian 
diffusion  process  described  by  the  ltd  stochastic  equation 

dx(t)  =  F_x(t)  +  Z1//2  dwt  . 

In  this  equation,  w  is  a  standard  N-dimensional  Wiener  process, 

F  =  -|j-  and  I  =  cov(£|f_(k)  =  f*).  (2.4) 

★ 

f 

where  £.  =  G.(f  .(k),  Y.(k),  U.(k)).  The  invariant  measure  for  (_x(t)) 

J  J  J  J  J 

is  Gaussian  with  mean  zero  and  covariance  matrix  Q  satisfying 

QF  +  FQ  +  I  =  0  .  (2.5) 

If  the  diffusion  approximation  is  sufficiently  accurate,  the 
covariance  matrix  of  f(k)  for  large  k  should  thus  be  Qe  +  o(c). 

Information  on  the  asymptotic  convergence  rate  of  solutions  to 
the  ODE  is  given  by  the  eigenvalues  of  F.  The  more  negative  the  real 
parts  of  the  eigenvalues,  the  quicker  the  convergence  rate.  The 
matrix  Q  gives  a  measure  of  steady  state  accuracy  of  the  control 
strategy.  As  a  specific  case,  (a"^Qa)e  gives  the  approximate  steady 
state  variance  of  a*f,  for  any  constant  N  vector  a. 
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Strategies  are  considered  such  that 

Gjd)  =  (fj)r  l^(f)  -aA^fj))  0<f  (2.6) 

where  r  and  a  are  positive  constants  with  r  >  1  and  a  >  o,  and  the 
following  conditions  are  satisfied  for  all  j: 

(condition  1)  ,$^(_f)  is  strictly  decreasing  in  f  and  is 
continuously  differentiable. 

(condition  2.j)  A.(0)  =  0,  A.(u)  is  strictly  increasing 

J  J 

in  u  and  is  continuously  differentiable, 
and  urA.(u)  has  a  finite  right-hand  derivative 

J 

at  u  =  0. 

These  conditions  ensure  that  G  is  locally  Lipschitz  continuous  so 
that  the  ODE  (2.3)  has  a  unique  solution  for  any  initial  point.  The 
restricted  region  of  ODE  (2.3)  considered  is 

-  =  (f  :  0  <  f.  <  A.  for  each  j) 

J  J 

Since  f.  ^cf-  +  o(f.)  for  all  j  and  for  some  negative  c,  the  solution 
J  J  J 

f(t)  remains  in  2  for  all  finite  t.  Therefore,  the  auxiliary 
variables  in  (2.3)  are  identically  zero. 

In  [3],  a  Liapunov-style  proof  of  global  stability  for  the  ODE 
(2.3)  is  presented.  This  proof,  however,  is  omitted  from  this  paper. 

Two  control  strategies  are  presented  which  satisfy  the  conditions 
set  forth  earlier  in  this  section.  A  detailed  explanation  of  these 
strategies  is  given  in  [3],  so  all  that  is  given  here  is  a  brief 
explanation  of  the  two  strategies.  The  first  strategy  considered  uses 
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the  common  O-l-e  feedback  to  adjust  the  transmission  probabilities. 

For  0-1-e  feedback,  Y.(k)=Z(k)  for  each  j,  where  Z(k)  is  equal  to  0,  1 

J 

or  e  depending  on  whether  zero,  one,  or  more  than  one  station  transmits 
in  slot  k.  The  transmission  rule  is  of  the  form  in  Eq.  (2.1)  where 


Gj(f,  y,  u)  =  (fj)[c(y)  -  oA^f^)] 


(2-7) 


for  the  vector  c_  =  (.418,  0,  -.582)  which  is  the  vector  suggested  in 

[3]  on  the  basis  of  a  Poisson  approximation.  The  arbitrary  functions 

A.  also  must  satisfy  condition  2.j  for  all  j. 

J 

The  second  strategy  considered  uses  the  acknowledgment  based 
transmission  control  to  determine  the  transmission  probabilities. 

This  strategy  relies  only  on  the  knowledge  of  whether  or  not  each 
station's  own  transmissions  are  successful.  The  update  rule  given  in 
[3]  is  (using  Z(k)  as  above) 


Vk+D  * 

where  A  =  1  - 


c  •  (u,Z,U) 


[fjOO  +  of  j  (k)  c  j  ( f  j  ( k ) ,  Z(k),  IL'k))]- 


(l-u)e  -  1  -  ,A.(u)  Z=l,  U=1 


■1  -  jA.(u) 


Z=e,  U=1 


and  the  functions  A.  are  arbitrary  functions  satisfying  condition  2.j 

J 

for  all  j.  Both  strategies  given  here  are  simulate!  under  th~ 
conditions  given  in  section  (2.A.I.),  and  the  results  of  tneso 
simulations  are  Given  in  the  final  two  sections. 
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The  results  of  the  simulations  are  compared  to  predictions 

determined  as  follows.  To  simplify  the  simulations,  all  of  the 

functions  Ai  are  the  same  for  all  i.  This  results  in  the  globally 

★ 

stable  equilibrium  point  f  of  the  ODE  (2.3)  having  the  form 
★ 

f  =  (u,  ...  ,  u),  where  u  is  easily  computed.  The  F  and  I  matrices 
defined  in  (2.4)  have  the  form 

F  =  ai  +  bW,  I  =  cl  +  dW  (2.8) 

where  a,  b,  c,  d  are  readily  computed  and  W  is  the  matrix  of  all  ones, 

and  so  the  scaled  covariance  matrix  Q  determined  by  (2.5)  also  has  the 

form  Q  =  rl  +  sW,  where  r  and  s  can  be  readily  expressed  in  terms  of 

a,  b,  c,  and  d.  In  general,  a  and  b  are  less  than  or  equal  to  zero. 

F  has  eigenvalue  a+b  with  multiplicity  one  (corresponding  to  the 

eigenvector  of  all  ones),  and  the  eigenvalue  a  with  multiplicity  (n- 

1).  The  smallest  magnitude  eigenvalue  of  F  is  used  as  an  indicator  of 

how  quickly  the  individual  transmission  probabilities  should  fall  into 

line,  or  of  the  relaxation  time  seen  in  the  system.  The  predictions 

given  in  [3]  for  the  "relaxation  time"  as  well  as  the  expected 

★ 

throughput  and  the  equilibrium  value  of  _f  are  compared  to  simulation 
values  in  the  final  two  sections  of  this  paper. 

The  control  of  the  sum  of  transmission  probabilities  is  also  used 
to  determine  how  well  simulation  values  conform  to  predictions.  The 
sum  is  written  as  e*f  where  e  is  the  all  ones  row  vector.  An  estimate 
for  the  steady  state  variance  of  each  f^  is  found  by  locking  at  the 
diagonal  entries  of  the  matrix  Q  and  multiplying  this  by  ~.  e‘Qe 
aporoximates  I/-:  times  the  steady  state  variance  of  e-f.  Also,  the 
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magnitude  of  the  eigenvalue  of  the  matrix  F  for  eigenvector  e  gives  a 
speed  of  convergence  measure  for  _e*f\  These  measures  are  used  to 
determine  if  the  simulated  models  follow  the  predicted  values  found 
using  the  theory  given. 

2.B.  Simulation  of  Nonzero  e,  Fixed  N 

2.B.I.  Models  and  Simulation 

The  models  used  in  this  section  are  specific  examples  of  the 
models  given  in  the  previous  section.  In  this  section,  N  is  fixed; 
however,  e  is  not  allowed  to  tend  to  zero  as  is  the  case  in  the 
previous  section.  The  update  rules  are  the  same  as  those  given  in  the 
previous  section  with  A^(u)=u  for  all  j.  The  choice  of  the  *  and 
values  is  discussed  in  the  following  subsections  along  with  the 
simulation.  These  simulation  results  are  then  compared  to  the 
predictions  made  in  [3]. 

A  Fortran  program  is  used  to  simulate  the  strategies.  The 
simulations  use  a  uniform  random  number  generator  to  determine  which, 
if  any,  stations  transmit.  Then,  transmission  probabilities  are 
updated  according  to  the  strategy  used  and  the  feedback  received.  The 
program  prints  out  the  throughput  and  the  transmission  probabilities 
for  each  station  at  predetermi ned  intervals  throughout  the  simulation. 
An  estimate  for  the  average  transmission  probability  is  determined 
over  the  total  number  of  slots  excluding  the  first  thousand  slots  in 
order  to  reduce  transient  effects.  Various  initial  conditions  are 
used  ror  the  simulations.  The  initial  transmission  probabi 1 i ties  are 
chosen  randomly  for  the  most  part.  In  tr.is  way,  predictions 
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concerning  the  transient  period,  such  as  the  relaxation  times,  can 
be  checked.  An  estimate  for  the  variance  of  the  individual 
transmission  probabilities  and  of  the  sum  of  transmission 
probabilities  using  the  same  sample  space  as  that  used  to  determine 
the  mean  is  also  evaluated.  All  of  these  estimates  are  printed  out  at 
the  end  of  a  simulation  run  and,  for  some  of  the  runs,  discussed  in 
the  next  subsection. 

2.B.2.  Results  of  the  N  Fixed  Case 

The  predictions,  from  [3]  and  the  estimates  from  the 
simulations  are  compared  in  this  subsection.  This  comparison  is  used 
to  determine  if  the  strategies  perform  as  expected,  and  hence  provide 
proof  that  the  simulations  bear  out  the  theory.  All  simulations  are 
run  over  4000  to  5000  slots,  the  only  exceptions  being  in  cases  noted 
in  this  subsection  where  the  number  of  slots  was  allowed  to  be  10,000 
and  30,000  slots. 

The  simulations  for  the  0-1-e  strategy  all  exhibit  the  fact  that 
once  two  or  more  stations'  transmission  probabilities  become 
approximately  equal  (i.e.,  fall  into  line),  the  probabilities  remain 
equal  throughout  the  remainder  of  the  simulation.  The  transmission 
probabilities  are  adjusted  by  the  same  update  rules  since  the  same 
feedback  is  used  by  each  of  the  stations.  This  is  unlike  the 
acknowledgment  based  strategy  where  each  of  the  stations  updates  its 
own  transmission  probabi 1 i ties  using  its  own  feedback  determined  by 
whether  or  not  it  transmits  in  a  slot. 

Table  (2.1)  shows  the  throughout  and  Table  (2.2)  shows  the  estima 
mean  transmission  probabi  1  i  ties  for  various  values  of  i  and  11,  and  a 
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Table  (2.1)  Throughput  results  for  the  0-1-e  strategy 
for  c=.05  and  for  various  N  and  j.,  given 
in  successful  transmissions  per  slot 


N=2 

N=6 

N=10 

=  .088 

expected 

simulated 

.49 

.47 

.40 

.39 

— 

=  .25 

expected 

simulated 

.46 

.47 

.40 

.39 

_  _  _  _ 

II 

o 

expected 

simulated 

.39 

.39 

.39 

.40 

.38 

.38 

=  1.0 

expected 

simulated 

.35 

.35 

.38 

.37 

.38 

.37 

Table  (2.2)  Transmission  probabilities  for  the  0-1-e 


strategy  for 

s=.05  and 

N=2 

N=6 

088 

expected 

.42 

.16 

simulated 

.41 

.14 

25 

expected 

.37 

.15 

simulated 

.36 

.15 

707 

expected 

.26 

.13 

simulated 

.26 

.13 

.0 

expected 

.22 

.  12 

simulated 

.22 

.11 
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value  of  .05  for  e.  These  values  are  compared  to  predicted  values  in 
the  same  tables.  In  all  cases  considered,  the  simulated  values  are 
near  the  predicted  values.  The  error  between  the  simulated  and  the 
predicted  values  is  given  by  the  steady  state  accuracy.  This  steady 
state  accuracy  is  measured  by  the  variance  about  the  mean  of  the 
transmission  probabilities.  This  is  discussed  in  more  detail  later  in 
this  subsection.  The  lack  of  results  for  the  N= 1 0  cases  where  a=.088 
and  a=.25  is  due  to  the  convergence  problems  present  when  a  small  a  is 
chosen.  As  a  gets  smaller,  the  smallest  magnitude  eigenvalue 
decreases  in  magnitude  for  a  constant  N.  This  results  in  a  longer 
relaxation  time  which  requires  a  longer  period  of  time  prior  to  the 
recording  of  data  due  to  a  longer  transient  period.  This  in  turn 
requires  more  memory  and  CPU  time  from  the  computer  which  soon  becomes 
expensive.  For  this  reason,  and  the  fact  that  the  simulations  done 
produce  adequate  proof  that  the  theory  predicts  the  simulation  results 
for  these  cases,  the  N  =  1 0  runs  are  omitted. 

We  now  consider  the  effect  of  changing  the  values  of  a  and  of  c 
on  the  steady  state  variance  and  the  relaxation  time.  As  a  and  i.  are 
increased  it  should  take  less  time  to  converge  to  the  predicted 
transmission  probabilities.  However,  this  causes  the  variance  of  the 
transmission  probabi I i ties  to  increase  which  is  equivalent  to  saying 
the  steady  state  accuracy  is  degraded.  On  the  other  hand,  an  and  an 
that  are  smaller  in  magnitude  result  in  a  longer  relaxation  time  which 
in  turn  causes  a  longer  transient  period  prior  to  convergence  as  in 
the  N= 1 0  case  discussed  earlier  in  this  subsection.  The  effects  of 
varying  the  __  and  e_  va'ues  are  much  more  dramatic  in  the 
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acknowledgment  based  strategy  where  the  transmission  probabilities 
never  actually  become  or  stay  equal  over  any  length  of  time,  resulting 
in  larger  variances  and  in  larger  relaxation  times. 

This  leads  to  the  results  presented  in  Tables  (2.3),  (2.4),  and 
(2.5).  Table  (2.3)  produces  a  comparison  of  the  predicted  variances 
and  the  estimated  variances  found  in  the  simulations.  Again,  in  most 
cases  the  estimates  from  the  simulations  compare  well  with  the 
predicted  values.  Table  (2.4)  presents  predicted  relaxation  times 
found  using  the  theory  presented  in  [3].  Note  that  for  N=10,  the 
relaxation  times  for  a=.088  are  approximately  2400  slots.  For  this 
reason,  the  a=.088  and  N=10  case  is  omitted.  The  N=6  and  a=.088  case 
also  has  a  long  relaxation  time;  however,  a  10,000  slot  run  is  done 
for  this  case.  Only  one  10,000  slot  run  is  performed  due  to  the  cost 
of  producing  such  simulations.  It  is  also  noted  that  for  N=2,  the 
relaxation  times  all  lie  below  600  slots.  In  the  runs  done  for  N=2, 
for  the  most  part,  the  time  needed  for  the  transmission  probabilities 
to  converge  to  a  common  value  is  well  within  the  600  slot  relaxation 
times  given.  The  relaxation  times  are  also  dependent  upon  the  initial 
transmission  probabilities  used  by  the  stations.  However,  no  analysis 
is  done  on  this  transient  period  to  determine  the  effects  of  various 
initial  conditions.  Rather,  the  only  real  interest  in  this  transient 
period  concerns  ensuring  that  the  transmission  probabilities  do 
converge  under  random  initial  conditions  to  the  predicted  values, 
which  in  all  cases  they  do.  In  Table  (2.5),  the  estimated  variances 
found  for  various  ;  are  given  with  N=6  and  0.  Using  the  theory 
given  earlier,  it  is  known  that  the  variances  and  the  relaxation  times 
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Table  (2.3)  Steady  state  variance  of  individual  transmission 
probabilities  for  the  0-1-e  strategy  where  £=.05 


and 

for  various  N 

and  a 

N=2 

N=6 

N  =  1 0 

p 

II 

o 

co 

CO 

expected 

.001 

.0002 

_ 

simulated 

.001 

.0005 

— 

a=.  25 

expected 

.001 

.0002 

— 

simulated 

.001 

.0002 

— 

P 

II 

o 

expected 

.0003 

.0001 

.0001 

simulated 

.0003 

.0001 

.0001 

a-1.0 

expected 

.0002 

.0001 

.0001 

simulated 

.0002 

.0001 

.0001 

Table  (2.4)  Relaxation  times  for  the  0-1-e  strategy 
where  e=.05  and  for  various  N  and  « 


N=2 

N=6 

N=10 

a=. 088 

540 

1500 

2400 

a=.  25 

220 

550 

880 

a=. 707 

110 

230 

340 

ct=l .  0 

90 

170 

260 

Table  (2.5) 

Variance  for 

various  e  and  N=6, 

ci=  1.0 

£-.5 

c=.05 

£=.005 

expected 

.001 

.0001 

.00001 

s imul a  ted 

.001 

.0001 

.00001 

for  that  matter,  are  directly  affected  by  the  value  chosen  for  e.  The 
choice  of  £=0.05  for  the  runs  already  discussed  is  a  value  originally 
chosen  after  some  experimentation  with  various  values  of  e.  This 
value  proved  to  be  sufficient  for  the  simulations  done,  but  by  no 
means  is  it  assumed  to  be  the  best  choice  of  e  possible. 

The  final  results  for  the  0-1-e  strategy  are  summed  up  in  Tables 
(2.6)  and  (2.7).  In  Table  (2.6),  the  predicted  and  estimated  variance 
of  the  sum  of  the  transmission  probabilities  is  given.  This  measure 
of  steady  state  accuracy  provides  a  means  to  see  how  well  as  a  whole 
the  transmission  probabilities  have  converged  to  the  predicted  mean 
sum.  Again,  the  results  appear  favorable.  Table  (2.7)  provides  the 
relaxation  times  for  the  sum  of  the  transmission  probabilities.  This 
presents  an  idea  as  to  how  quickly  the  sum  converges  for  the  parameter 
values  discussed. 

This  leads  to  the  discussion  of  results  for  the  simulation  of  the 
Ack  strategy.  Tables  (2.8),  (2.9),  (2.10),  (2.11),  (2.12),  and  (2.13) 
present  the  information  analogous  to  that  presented  in  the  previous 
tables  for  the  0-1-e  strategy.  In  Table  (2.8),  the  predicted  throughput 
and  estimated  throughput  are  given.  Again,  the  estimates  follow  the 
predictions  for  each  of  the  cases.  Table  (2.9)  presents  the  predicted 
mean  transmission  probability  and  the  simulated  estimated  mean  with  a 
range  of  values  the  individual  means  took  at  the  conclusion  of  the 
simulation.  For  the  Ack  strategy,  as  discussed  earlier,  the 
individual  stations  have  transmission  probabilities  which  do  not 
ultimately  stay  equal.  For  this  reason,  the  estimated  means  of  the 
individual  stations  take  on  a  range  of  values.  The  theory  presented 


Table  (2.6)  Variance  of  sum  for  0-1-e  strategy  for 
£=.05  and  various  N  and  a 


N=2 

N=6 

N=10 

a= . 088 

expected 

.005 

.008 

simulated 

.004 

.008 

— 

a=.  25 

expected 

.003 

.007 

— 

simulated 

.005 

.007 

---- 

a=. 707 

expected 

.001 

.004 

.006 

simulated 

.001 

.005 

.006 

o=1.0 

expected 

.001 

.004 

.005 

simulated 

.001 

.003 

.005 

Table  (2.7)  Relaxation  times  of  the  sum  of  transmission 
probabilities  for  e=.05  and  various  N  and  a 


N=2 

N=6 

N*10 

P 

II 

o 

00 

CO 

80 

260 

430 

a=.  25 

80 

210 

340 

a=. 707 

70 

140 

220 

[•; 


Table  (2.8)  Throughput  for  Ack  strategy  for  £=.05  and 
various  values  of  N  and  a,  given  in 
successful  transmissions  per  slot 


«=  088 

expected 

N=2 

.47 

N=6 

.40 

N=10 

•  WOO 

simulated 

.50 

.44 

— 

a=.  25 

expected 

.46 

.40 

—  —  •  — 

simulated 

.47 

.40 

— 

P 

It 

O 

^sl 

expected 

.44 

.40 

.39 

simulated 

.44 

.39 

.39 

p 

II 

h- * 

o 

expected 

.43 

.38 

.38 

simulated 

.40 

.38 

.38 

Table  (2.9)  Transmission  probabilities  for  Adk  strategy  for 
£=.05  and  various  values  of  N  and  a  (with  an 
average  of  estimates  and  a  range  of  estimates 
for  the  simulated  case) 


N=2 

N=6 

N  =  10 

« 

ii 

o 

cx> 

00 

expected 

simulated 

.38 

. 35( . 21- . 50) 

.15 

.  13( . 03- . 26) 

— 

a=.  25 

expected 

simulated 

.37 

. 36( . 33- . 39) 

.15 

.  14( .06- . 20) 

_  _  _  _ 

P 

II 

'vl 

O 

expected 

simulated 

.33 

. 32( . 30-  .  33) 

.14 

. 13( . 08- .19) 

.09 

. 08( . 008- .17) 

3*1.0 

expected 

simulated 

.31 

. 29( . 25- . 33) 

.14 

.  12 ( . 07- .17) 

.09 

. 08( . 03-  .16) 

Table  (2.10) 


Steady  state  variance  of  individual 
transmission  probabilities  for  the  Ack 
strategy  where  e=.Q5  and  for  various  N  and  a. 
(Note  that  the  simulated  variance  is  an  average 
of  all  the  estimated  variance  and  a  range  of 
values  is  also  given. ) 

N=2  N=6  N=10 


a=.088 

expected 

simulated 

.02 

.005( .004-. 005) 

.04 

.008 (.00008-. 02) 

.... 

a=.  25 

expected 

simulated 

.006 

. 003{ .002- .004) 

.015 

.001 (.0002-. 003) 

a=. 707 

expected 

simulated 

.003 

•003( .002-. 003) 

.005 

.001 (.0005-. 003) 

.004 

.001(.000-.005) 

a=l .  0 

expected 

simulated 

.002 

•002(.002-.002) 

.004 

.002( .0002-. 004) 

.003 

.0006 (.00005-.  001) 

Table  (2.11) 


Relaxation  times  of  individual  transmission 
probabilities  for  the  Ack  strategy  where 
e=.05  and  for  various  N  and  a 


N=2 

N=6 

N=10 

p 

m 

o 

00 

00 

1500 

9900 

25000 

in 

C\J 

ii 

600 

3600 

9200 

ii 

o 

270 

1400 

3500 

a=l .  0 

210 

1100 

2600 

oviVs  kk  o-o.  ..n 


Table  (2.12)  Variance  of  sum  for  Ack  strategy  for  e=.05 
and  various  N  and  a 


N=2 

N=6 

N=10 

a=. 088 

expected 

.006 

.007 

_ 

simulated 

.004 

.040 

— 

a=.  25 

expected 

.006 

.006 

— 

simulated 

.006 

.006 

— 

P 

II 

o 

expected 

.005 

.006 

.004 

simulated 

.004 

.005 

.005 

p 

II 

►— * 

o 

expected 

simulated 

.004 

.004 

.005 

.005 

.004 

.004 

Table  (2.13) 

Relaxation  times  of 
probabilities  for  e= 

the 

.05 

sum 

and 

of  transmission 
various  N  and  a 

N=2 

N=6 

N=  1 0 

00 

00 

o 

II 

n 

80 

680 

1900 

ct=.  25 

80 

620 

1700 

a=. 707 

70 

510 

1300 

a=l .  0 

70 

470 

1200 
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earlier  predicts  this  result  as  the  data  in  Tables  (2.10)  and  (2.11) 
illustrate.  Table  (2.10)  contains  the  predicted  variances  and  the 
simulated  range  of  variances  for  the  various  cases.  It  is  seen  that 
these  variances  are  quite  large  compared  to  the  variances  of  the  0-1-e 
strategy.  Hence,  it  is  expected  that  the  various  transmission 
probabilities  vary  more  about  the  mean  than  in  the  0-1-e  strategy. 

Table  (2.11)  presents  the  predicted  relaxation  times  again  found  from 
the  theory  discussed  earlier.  It  is  noted  that  these  times  are  much 
greater  than  the  times  for  the  0-1-e  strategy.  Hence,  it  is  expected 
that  a  larger  transient  period  also  influences  the  results  already 
presented.  Since  all  the  runs  already  discussed  are  over  relatively 
short  simulation  times,  10,000  slots  or  less,  a  run  of  30,000  slots  is 
done  to  minimize  transient  effects,  al though  the  results  given  thus  far 
for  the  simulations  done  do  follow  the  predictions  from  the  theory 
consistently.  The  run  of  30,000  slots  for  N=10,  a=.707,  and  e=.Q5  is 

discussed  later  in  this  subsection. 

Tables  (2.12)  and  (2.13)  contain  the  analogous  data  to  that  given  fo 
the  variances  of  the  sum  and  relaxation  times  of  the  sum  for  the  0-1-e 
strategy.  Again,  the  simulation  results  are  consistent  with  the 
predictions.  It  is  worth  noting  that  these  results  imply  that  the  sum 
converges  rather  quickly,  compared  to  the  individual  station's 
relaxation  times,  to  a  certain  value;  however,  it  is  known  from  the 
given  results  for  the  individual  stations  that  the  individual 
transmission  probabilities  vary  to  a  much  larger  extent.  This  implies 
that  a  judgment  as  to  how  well  a  simulation  is  following  predictions 
cannot  be  made  by  looking  only  at  the  throughput  and  results 
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concerning  the  sum  of  transmission  probabi 1 ities.  So,  to  alleviate  any 
problems  over  ambiguity  of  results  and  over  terribly  long  transient 
periods,  a  30,000  slot  simulation  is  done. 

The  parameters  used  for  this  run  are  e=.05,  a=.707,  and  N=10. 

The  choice  of  30,000  slots  is  due  to  limitations  in  the  computer 
system  employed  to  do  this  simulation.  This  allows  ample  time  for 
recoveries  by  all  the  stations  since  the  predicted  relaxation  time  is 
3500  slots.  Table  (2.14)  contains  some  of  the  data  found  for  this 
simulation  and  compares  these  to  the  predicted  values.  Again,  it  is 
seen  that  the  simulation  results  are  consistent  with  the  predictions. 
This,  however,  says  very  little  about  whether  individual  stations 
recover  from  cases  where  a  station's  transmission  probability  falls 
well  above  or  well  below  the  expected  value.  To  say  something  about 
this  phenomenon,  a  definition  of  recovery  is  needed.  A  rough 
definition  of  recovery  used  in  this  paper  is  that  a  recovery  occurs 
when  a  station's  transmission  probability  returns  from  a  value  outside 
of  plus  or  minus  one  standard  deviation  about  the  mean  to  a  value 
within  one-half  the  standard  deviation  about  the  mean.  This 
definition  is  used  here  only  as  a  measure  of  recovery  for  this 
simulation.  Other  possible  definitions  of  recovery  exist;  however, 
this  is  the  definition  adopted  here. 

Table  (2.15)  presents  a  short  excerpt  from  the  data  obtained  from 
the  simulation.  These  are  just  five  examples  of  recoveries  from  a 
total  of  thirty-six  recoveries  which  occurred  using  the  definition 
discussed  earlier.  Again,  the  predictions  hold  true  first  in  allowing 
fo-'  stations  to  fall  outside  of  the  recovery  region  and  second  in 


Table  (2.14)  Results  for  the  N=10,  30,000  slot  simulation 


Expected 

Simula  ted 

Throughput 

.39 

.39 

Transmission  Probabilities 

.09 

.08  (.04-. 13) 

Variance 

.004 

.002  (.0003-. 003) 

Relaxation  time 

3500 

— 

Variance  of  sum 

.004 

.005 

Relaxation  time  of  sum 

1300 

_____ 

Table  (2.15)  Examples  of  recovery  times  for  various  stations 
under  the  definition  of  recovery  given  in  the 
text 

Transmission  Probabi 1 i ties 


Slot  Number 


Station  s 


Station  4 


Station  10 
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predicting  relaxation  times  that  most  of  the  recoveries  fell  within. 
The  average  recovery  time  over  the  thirty-six  recoveries  is  found  to 
be  3200  slots  for  the  simulation  done.  The  final  set  of  predictions 
concerns  those  made  for  the  variance  of  the  sum  and  the  relaxation 
time  of  the  sum.  These  results  are  also  presented  in  Table  (2.14). 
Again,  the  predictions  are  consistent  with  the  relaxation  time  found 
for  the  simulation  falling  within  the  predicted  value. 

All  these  results  supply  evidence  that  the  theory  developed  in 
[3]  and  presented  in  a  condensed  form  in  an  earlier  subsection 
provides  consistent  predictions  as  to  the  behavior  of  the  control 
strategies  with  nonzero  e.  Although  little  is  said  about  how  useful 
these  strategies  are  in  reality,  this  has  nothing  to  do  with  the  main 
result  that  it  is  possible  to  predict  the  performance  of  the 
strategies  with  the  given  theory.  The  simulations  provided  some  proof 
of  this.  As  a  further  variation  and  test  of  the  power  of  the  theory, 
predictions  concerning  the  performance  of  a  varying  N  case  are 
compared  to  simulation  results  as  done  in  this  subsection. 

2.C.  N  Allowed  to  Slowly  Vary 

2.C.I.  Models  and  Simulation 

The  model  discussed  in  this  subsection  differs  from  the 
model  presented  in  the  previous  section  only  in  that  this  model  allows 
the  number  of  stations,  N,  to  vary  in  the  simulation.  The  control 
strategies  and  the  basic  programming  format  are  the  same  as  those 
given  in  subsection  2.B.I. 
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Two  different  methods  are  employed  to  vary  N.  The  first  is  a 
deterministic  strategy  where  the  number  of  stations  is  given  by  a  sine 
function  as  follows: 

N(s)  =  A  +  B  Integer  [.1  +  sin(Cs)]  (2.9) 

where  N(s)  is  the  number  of  stations  used  during  slot  s;  A  and  B  are 
constants  determining  the  range  of  N(s);  C  is  a  constant  in  radians 
per  slot  denoting  the  frequency  of  change,  and  Integer  [  ]  takes  the 
integer  part  of  the  expression  in  the  brackets.  This  defines  the 
deterministic  system.  Although  the  number  of  stations  and  the  points 
in  time  where  the  number  of  stations  change  are  known,  the  choice  of 
users  to  be  used  for  a  particular  length  of  time  is  made  randomly. 

For  example,  consider  the  case  where  there  are  six  particular  stations 
out  of  a  possible  15  stations  using  the  channel  for  the  next  30  slots, 
and  after  that  the  number  of  stations  is  increased  to  seven  for  the 
next  20  slots.  The  seven  stations  need  not  include  all  the  six 
previous  stations  since  at  each  change  it  is  equally  probable  that 
each  of  the  individual  stations  will  be  used  in  the  next  set  of  time 
slots. 

The  second  method  employed  chooses  both  the  number  of  stations 
and  the  number  of  slots  these  stations  are  "on"  randomly  from  a  set  of 
possible  stations  and  slot  times.  This  method  is  called  the  random 
system.  The  choice  of  the  number  of  stations  used  is  dependent  on  a 
uniformly  distributed  random  variable  which  takes  on  values  from  0  to 
1.  The  function  used  to  determine  11  is  given  as  follows: 


N  -  K  +  Integer  [L  +  R] 


(2.10) 
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where  K  and  L  denote  constants  determining  the  range  of  N.  R  is  a 
uniformly  distributed  random  number  as  discussed  above.  Again  the 
choice  of  which  N  stations  that  are  "on"  for  a  certain  period  is  done 
randomly  from  the  set  of  total  possible  stations.  The  period  time,  or 
the  amount  of  time  prior  to  changing  the  number  of  stations  again,  is 
determined  by  another  uniformly  distributed  random  variable  which 
takes  on  values  from  0  to  1.  This  function  is  given  as  follows: 

s  =  D  +  Integer  [E  R]  (2.11) 

where  D  and  E  determine  the  range  of  s,  and  R  is  the  same  as  defined 
above. 

Another  choice  made  for  the  simulations  is  the  choice  of  initial 
conditions  for  each  station.  Each  station  transmits  with  probability 
1  for  the  0-1-e  strategy  and  (1-e*1)  for  the  Ack  strategy  when  it 
initially  turns  on.  However,  two  different  methods  are  used  when  the 
stations  are  turned  on  in  the  simulation  after  having  already  been 
turned  on  and  then  off.  The  transmission  probability  is  set  to  one 
for  the  0-1-e  strategy  or  2  (1-e’^)  for  the  acknowledgment  based 
strategy  in  the  first  method.  The  second  method  allows  each  station  to 
save  its  final  transmission  probability  prior  to  being  turned  off, 
allowing  it  to  be  used  again  when  that  station  is  turned  on  again. 
These  methods  are  used  only  to  gain  insight  into  how  to  choose  initial 
conditions  under  these  circumstances.  No  analysis  is  done  to 
determine  which  method  is  superior. 

When  the  number  of  stations  change,  the  throughout,  up  to  that 
point  is  printed  out.  The  number  of  stations  and  the  time  slot  are 


also  printed  out.  The  throughput  given  from  the  simulation  is  then 
compared  to  the  predicted  throughput  from  the  theory  in  [3]  to 
determine  how  well  the  simulation  is  predicted.  Throughput  may  not  be 
the  most  informative  measure  of  performance  since  the  individual 
transmission  probabilities  are  not  compared  to  their  predicted  values. 
The  throughput  is  a  simple  measure  that  tests  the  predictions  for  the 
case  that  N  is  slowly  varied. 


2.C.2.  Comparison  of  the  Predicted  Throughputs 
and  of  the  Simulated  Throughputs 

The  simulated  throughputs  are  summarized  for  both  strategies 
in  Table  (2.16)  for  the  specified  conditions.  These  throughputs  are 
found  using  the  values  given  for  the  parameters  in  equations  (2.9), 
(2.10),  and  (2.11)  as  follows: 


A  =  10 


B  =  5 


C  =  2  / 300 
D  =  20 


E  =  100 


U  =  5 


L  =  11 


and  using  600  slot  runs. 

The  theory  yields,  for  given  e>  ot,  and  average  N  values,  a 
relaxation  time  such  that  if  N  is  slowly  varying  in  a  time  scale  on 
the  order  of  the  relaxation  time  then  the  transmission  probabilities 
are  controlled  to  within  a  predicted  tolerance.  For  the  systems 
discussed,  an  average  N  of  10  is  used  as  the  N  needed  to  predict  the 
relaxation  times.  In  Table  (2.17),  the  relaxation  times  that  cannot 
be  exceeded  if  the  theory's  predictions  are  to  hold  are  presented. 
Also,  in  Table  (2.17)  there  are  the  predicted  throughputs  for  the 
cases  that  are  simulated. 


Table  (2.15) 


Table  (2.17) 


Throughputs  found  for  the  variable  N  simulations 
where  ct= .  25  and  the  runs  covered  6000  slots. 

Note:  _  has  reinitialized  to  one  policy  on  top 
/  #  has  rei nitial ized  to  old  value  policy 
on  bottom 

Random 


Ack 

0-1 -e 

.28 

' .  34 

.26, 

'.25 

.35, 

'.39 

.36, 

'.36 

Deterministic 

.30, 

'.35 

.33, 

'.28 

.39, 

'.34 

.37 

'.36 

Predicted  throughputs,  relaxation  times,  and 
relaxation  times  of  the  sum  for  N=10,  a  =  .25 

Relaxation  Time 

Throughput  Relaxation  Time  at  Sum 


The  random  system  is  designed  to  have  an  average  period  of  70 
slots  prior  to  a  change  in  the  number  of  stations.  The  deterministic 
system  is  designed  so  that  the  number  of  slots  it  takes  to  go  from  5 
to  15  stations  is  150.  Hence,  the  models  of  varying  N  which  we  used 
in  our  simulation  have  N  roughly  constant  over  time  periods  on  the 
order  of  70  slots  for  the  random  system  and  100  slots  for  the 
deterministic  system.  These  time  periods  exceed  the  theoretical 
relaxation  times  in  only  one  of  the  four  different  cases,  namely  the 
case  where  e=.5  and  a=.25  for  the  0-1-e.  Otherwise,  the  N  is  varying 
faster  than  the  theory  allows,  so  under  these  cases  the  series  need 
not  apply  to  predicting  the  throughput.  A  solution  to  this  is  to  slow 
down  the  rate  of  variation  of  N  allowing  the  resulting  relaxation 
times  to  be  larger  than  what  is  required  for  the  theory. 

In  the  cases  that  the  theory  does  apply,  it  is  seen  that  the 
results  are  fairly  accurate  in  that  the  simulated  throughputs 
approximate  the  predicted  throughputs.  When  the  theory  does  not 
apply,  as  in  the  other  three  cases,  the  strategies  appear  to  still 
achieve  good  performance  (i.e.,  high  throughput).  Here,  a  more 
"global"  theory  is  needed  to  predict  performance  when  the  local  theory 
described  in  [3]  does  not  apply. 

It  is  noted  that  any  comparison  of  the  two  policies  as  to  how  to 
set  the  transmission  probabilities  when  a  station  turns  on  after 
already  being  on  once  produces  no  useful  results.  Neither  policy 
appears  to  be  better  than  the  other  using  the  data  in  Table  (2.16)  as 
the  measure.  Hence,  no  conclusions  are  raised  about  this  point. 
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